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ABSTRACT 

We show that the equations of motion governing the evolution of a colhsionless grav- 
itating system of particles in an expanding universe can be cast in a form which is 
almost independent of the cosmological density parameter, Q, and the cosmological 
constant, A. The new equations are expressed in terms of a time variable t = InD, 
where D is the linear rate of growth of density fluctuations. The weak dependence on 
the density parameter is proportional to e = — 1 times the difference between the 

peculiar velocity (with respect to r) of particles and the gravity field (minus the gra- 
dient of the potential), or, before shell-crossing, times the sum of the density contrast 
and the velocity divergence. In a 1-dimensional collapse or expansion, the equations 
are fully independent of fl and A before shell-crossing. In the general case, the effect of 
this weak f2 dependence is to enhance the rate of evolution of density perturbations in 
dense regions. In a flat universe with A 7^ 0, this enhancement is less pronounced than 
in an open universe with A = and the same f2. Using the spherical collapse model, 
we find that the increase of the rms density fluctuations in a low universe relative to 
that in a flat universe with the same linear normalization is ^ O.Ole(ri) <S^ >, where 6 
is the density field in the fiat universe. The equations predict that the smooth average 
velocity field scales like il"-^ while the local velocity dispersion {rms value) scales, ap- 
proximately, like ri'^'^. High resolution N-body simulations confirm these results and 
show that density fields, when smoothed on scales slightly larger than clusters, are 
insensitive to the cosmological model. Halos in an open model simulation are more 
concentrated than halos of the same M/fl in a flat model simulation. 

Key words: cosmology: theory - dark matter - large scale structure of Universe 



1 INTRODUCTION 

The cosmological background determines the growth rate of 
matter density fluctuations. This is the result of two effects. 
First, the initial conditions are specified in terms of the den- 
sity contrast field S = p{x)/pb — 1. Therefore the actual den- 
sity, p{x), which dictates the dynamical evolution, as can be 
seen for example from the spherical top-hat model, depends 
on the mean matter density, pt- The second effect comes 
about simply because the mean matter density varies with 
time according to the assumed cosmological model. This in 
turn translates into a dependence of the evolution of the fluc- 
tuation field S on the parameters of the cosmological model: 
the density parameter, SI, and the cosmological constant, 
A. Here we focus on the following aspect of the dependence 
of dynamics on the cosmological background. Starting from 
an initial density fluctuation field and a given amplitude 
of the evolved field, we address the question: how do the 
evolved peculiar velocity and density fields depend on the 
parameters and A? In the linear (e.g. Peebles 1980) and 
in the Zel'dovich quasilinear (Zel'dovich 1970) approxima- 



tions, once an initial density fluctuation fleld is evolved to 
a given amplitude, it does not contain any information on 
the parameters Q and A. In these approximations, the pe- 
culiar velocity fleld is simply proportional to f{Q,A) where 
/ is the so-called linear growth factor. This result is easy to 
understand. In the linear approximation, the density fluc- 
tuations are merely amplifled by a time dependent factor, 
D. In the Zel'dovich approximation, the displacement vec- 
tor is the product of the initial gravity fleld and the func- 
tion D. Moreover, second order perturbation theory calcu- 
lations (e.g. Bouchet et al. 1992) have shown that moments 
of the density fluctuation field are very insensitive to Q, and 
A. Finally, in the highly nonlinear regime, N-body simu- 
lations (e.g. Davis et al. 1985) show that the final matter 
distribution in simulations with the same initial conditions 
changes very little as the parameters of the cosmological 
background are varied. Significant differences between flat 
and open models are found only in the cores of what are 
identified as rich clusters in these simulations. These results 
have proved useful in analyzing observations of the large 
scale structure. Nusser and Dekel (1993), for example, used 
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N-body simulations to argue that a recovery of the initial 
density fluctuations from the observed galaxy distribution 
is almost independent of /(fi, A), where as a recovery from 
the observed peculiar velocity field is sensitive to the as- 
sumed Q,. They applied their reconstruction method to the 
POTENT compilation of the peculiar velocity data and to 
the 1.2 Jy IRAS survey and concluded that f2 > 0.3 with 
high confidence. Bernardeau et al. (1995) used second or- 
der perturbation theory to argue that the reduced skewness 
of the divergence of the peculiar velocity field is inversely 
proportional to /, in accordance with the scaling implied 
by the Zel'dovich approximation. They found that the mea- 
sured skewness is consistent with Q, of about unity. 

Here we aim at a better understanding of the depen- 
dence of the equations of motion on the cosmological pa- 
rameters. In section 2 we write the equations of motion in a 
form which is almost independent of the background cosmol- 
ogy. We discuss the Q. dependence in toy-models in section 
3. In section 4, we use high resolution N-body simulations to 
investigate in detail the differences in the matter distribu- 
tion in fiat and open models. We conclude with a summary 
and discussion in section 5. 



ALMOST Q. INDEPENDENT EQUATIONS 
OF MOTION 



yields Q. = cqIH^e? where co = Q,Qii^a% and the subscript 
denotes quantities at the present time. Therefore, energy 
conservation and (0) yield 



Co 



Co ~f -I -f fca 



(5) 



The bubble factor, i/, can be eliminated from these 
equations by working with a new time variable p = Ina. We 
define a new "velocity" a. = da;/dp = H~^v. In these new 
variables, the Poisson equation remains unaltered while the 
continuity and Euler equations become 



+ {1 + S)V ■a = 0, 



dS 
dp 

g + (l-,)a=-^^7V^, 



(6) 



(7) 



where q{p) — Q/2 ~ \ is the time dependent deceleration 
parameter and we have used d_ff/dp — —(1 + q)H to derive 
(@). 

Attempting to eliminate f2 from ^ we further make 
an additional transformation from the time variable p to r 
defined by 



T = lnD{p), 



(8) 



We restrict our treatment to the case of a matter dominated 
universe with a cosmological constant, i.e, we assume that 
the total mean density is ptot = pb + A/3 where pb{t) is 
the matter contribution and A is the cosmological constant. 
We use the standard notation in which, a{t) is the scale 
factor, H{t) = a/a is the time dependent Hubble factor, 
n = ph{t)/pc{t) and A = A/3H^ where pc = 3H^/8tvG is 
the critical density. Let x and V = dx/dt be the position 
and peculiar velocity of a particle in comoving coordinates. 
The equations of motion are: the continuity equation 



dS_ 

dt 



+ {1+S)V ■v = 0, 



the Euler equation of motion, 



(1) 



(2) 



where D is the linear growing density mode determined by 
the equation. 



dt^ dt 2 



(9) 



Analytic solutions to Q can be found in Heath (1977). Ex- 
pressing (^ in terms of r defining the "velocity" j9 = dx/dr, 
continuity equation is 



d5_ 
d7 



(l + <5)6' = 0. 



and the Euler equation is 



(10) 



(11) 



and the Poisson equation, 



= S. 



(3) 



Note that we have defined = 2^g/3Q where $3 is the pe- 
culiar gravity potential in comoving coordinates. Equations 
(0), (^ and (j^) together with the Friedman equations for 
the background quantities Q and a fully specify the dynam- 
ics of pressure free density fiuctuations. The scale factor, a, 
can be solved for using the Friedman equation 



day 8n 2 A 2 
-j ^-Gp,a +-a 



(4) 



where k = +1, —1 and correspond, respectively, to closed, 
open or flat universes. Energy conservation, pta^ — const, 



where, — — V-j9, and / = dr/dp is the linear growth factor 
which relates density contrast, S, to the divergence of the pe- 
culiar velocity field, v, in the linear regime. For A = 0, a good 
approximation for / is / « fi" *^ (Peebles 1980). For A 7^ 0, 
Lahav et al. (1991) found that / ^ fi"'^ A(l n/2)/70. 
Therefore, for reasonable values of the cosmological con- 
stant we neglect the dependence of / on A in the approxi- 
mate forms for / (Lahav et al. 1991). Note that the velocity 
— v/{Hf). In the Zel'dovich approximation, this is the 



'(!-§)/- 



* For A = 0, the function / satisfies (1 — 

= 0, which to first order in 1 — C yields / ^ H^/^ for C 1 (see 
also Lightman & Schechter 1990). However, the general solution 
to this equation is better fitted by / fs fZ" '' for Q < 0.7. A fit 

which works well for 0.05 < < 1 is / = Jn . 
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displacement vector of a particle from its initial to present 
position. Using (|^, we find 

%^i^-'i)f=l^^f^ (12) 
so the Euler equation ( |ll| ) is 

^--9-^[l + ^W](fl-^)=0, (13) 
where 

e{^) = - 1 ^ - 1. (14) 

and we have defined g = — V<^. The weak dependence on Q, in 
( p^ through couples to the difference between the ve- 
locity, /3, and the gravity field, g. Since 57 ~ f at early times, 
initially the function e almost vanishes, thus, any changes in 
the dynamics as a result of this weak dependence on Q. occur 
at later times. 

In virialised regions, the acceleration of a particle is 
dominated by the gravity field g. It is easy to see that by ne- 
glecting the terms involving the velocity in ( p^ and working 
with a new time variable with respect to which the velocity 
is (1 -I- e)~^^^fi « we obtain an Euler equation which 

is independent of Q.. This velocity is approximately equal to 
the comoving peculiar velocity divided by HOP'^ . This scal- 
ing with Q. is not surprising since the virial theorem implies 
that the rms value of the physical velocities in virialised re- 
gions with a given density contrast has similar scaling with 
Q.. Note however that the scaling is only approximate since 
the density profile of virialised objects has some dependence 

on n. 



Figure 1. The quantities <5 and d versus tiie time t for a positive 
perturbation for various values of and Ao as indicated in the 
figure. The upper steeply rising and the lower curves, respectively, 
correspond to bound and unbound perturbations in the Ho = 
0.2 cases. The values of 5 at the turnaround radii of the bound 
perturbations are 4.5, 6.7 and 11.5 for (Qq, Ao)=(l., 0), (0.2,0.8) 
and (0.2,0) respectively. 



3 O DEPENDENCE IN TOY MODELS 

It is clear from the form of (|l^) that the source term which 
drives the evolution is larger for lower Q.. Therefore we ex- 
pect to see more evolved clustering in a low 57 universe than 
in an n = 1 universe with the same initial conditions and 
linear normalisation. It is instructive to investigate the ef- 
fect of the term e in cases of special symmetry. Consider first 
the spherical expansion or collapse before the occurrence of 
shell crossing. This case illustrates the effect of of changing 
Q. in the before shell-crossing. In cases of special symmetry, 
we find it easier to solve directly for 5 and 6 rather than for 
g and j3. Therefore we take the divergence of (^^ and use 
the Poisson equation to obtain 

^-e-tf -|[l + e(n)](5-e)=0, (15) 
dr 2 

where = ji^^iPo f = ^V^/" with J\f = 1, 2 and 3 
at the centers of configurations with planar, cylindrical and 
spherical symmetry, respectivel y. T herefore, in the spherical 
top-hat model, the equations (|l^), (|l^) together with the 
equations relating a, r, 57 and / are sufficient to determine 



Figure 2. The same as figure 1 but for negative perturbations. 
The 5 curves are almost indistinguishable. 

the evolution of the quantities 9 and 5. For A = 0, the spher- 
ical collapse model can be solved analytically (e.g. Peebles 
1980) if the initial peculiar velocity is neglected. Here we 
numerically integrate the equations (14) and (17) under the 
initial conditions, 5i — 9i with \6i\ <^ 1, where the subscript 
i refers to quantities at the initial time. These initial con- 
ditions are realistic as they arise naturally in linear theory 
(e.g. Peebles 1980). 

In figures 1 and 2 we show the density contrast ver- 
sus the time r for positive {5i > 0) and negative {5% < 
0) tophat perturbations for three background cosmologies: 
(57o, Ao)=(l,0), (0.2,0.8) and (0.2,0). Figure 1 shows curves 
for two values of the initial density contrast corresponding to 
bound and unbound perturbations for the 57o — 0.2 cosmolo- 
gies. Although the equations were integrated from r = —5.4 
to 0, for the sake of clarity, figure 1 shows results only for 
r > — 1. For bound perturbations, the growth of 5 and 9 is 
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Figure 3. The density parameter versus the time t for an open 
universe (dashed line) and a fiat universe with a cosmological 
constant (solid line). 



fastest in the open universe case, (fio = 0.2, Ao =0). How- 
ever, significant deviations appear only when 5 is larger than 
10 or so. This is consistent with the work of Peacock and 
Dodds (1996) who found that nonlinear effects in the evo- 
lution of power spectra in N-body simulations are stronger 
in an open universe than in a flat -I- A universe of the 
same SI. Although the -I- A = 1 case shows more rapid 
evolution than the Q — 1 case, the corresponding curves are 
very similar even when the densities are larger than their 
values at turnaround. The reason for this is clear from fig- 
ure 3 which shows SI versus r. For Aq 7^ we see that Q is 
almost unity until relatively late times. Therefore, until late 
times, the evolution of S and 6 is very similar to the f2o = 1 
case. For positive unbound and negative perturbations, the 
effect of the cosmology on the evolution of S is almost negli- 
gible. The 6 curves show some differences. In voids, 6 grows 
more slowly in the low fio models once the density contrast 
approaches —1. Since perturbations in flat universes with a 
cosmological constant evolve similarly to those with A = 0, 
we do not discuss the case A 7^ further. 

The top-hat model can be used to evaluate how the vari- 
ance of an evolved generic density fluctuation field depends 
on n. We require here that the density field is smoothed 
on large enough scales such that shell-crossing is removed. 
Curves of 5(0, < 1) versus 5i = 5{Q = 1) satisfy the follow- 
ing (empirical) relation 



S{fl) = 5i exp 



Si 



(16) 



free of Q. Therefore, we replace A, in (|l|), with 4A/{Af~lf 
where A/" is the "dimensionality" of the collapse at each point 
in space. Other than for purely symmetric configurations, 
the quantity A/" is somewhat ambiguous. One possibility is 
to define it at any point in space in terms of the eigenvalues 
of the initial velocity deformation tensor, say, as the ratio 
of the square of the sum of the eigenvalues to the sum of 
their squares. Nevertheless, for our purposes it is not crucial 
to specify the form of Af and we simply treat it as a fac- 
tor which depends on the local topology of the density field. 
If P is the probability distribution function of the field Si 
evolved with = 1, then it can be shown that the variance, 
(j^(Sl), of the field S{Q) is given by 



\n) = / p{Si) 



l + Si 

[ + s{n) 



[1 + <5(!^)]^ d5i - 1, 



(18) 



where we have introduced the ratio {1 + 5i)/{l + S) to ac- 
count for the change in the volume element as a result of 
the change in density. By expanding ( p^ to third order in 
Si and substituting the result in (hsl) we find 



= i-f 



ajv A 



+ S3 



(i+a^A) 



.A2' 



(19) 



where is the variance of the field 5i and ajv = 4/ (A/"— 1)^ 
where < TV < 3 is some number describing the dimension- 
ality of typical collapse configuration. Note that, in deriving 
the last relation, we have neglected any local correlation be- 
tween A/" and 5i. The factor should depend on the power 
spectrum of the initial fluctuations and it can be determined 
empirically from N-body simulations. 

Consider now the effect of changing flo on the dynamics 
in shell-crossing regions in the case of 1-dimensional collapse. 
This case is particularly instructive since solutions {13 = g) 
to the 1-dimensional equations of motion are fully indepen- 
dent of until the occurrence of shell-crossing. Unfortu- 
nately, even in the simple one dimensional collapse, we have 
no analytic solutions in shell-crossing regions. Therefore we 
first use the Zel'dovich solution until the formation of the 
first singularity. Then we switch on to a one-dimensional 
N-body code to move particles according to ( [l^ ) in the 
shell- crossing phase. The initial density field we choose is 
Si oc cos(a;). Results of simulations with f2o = 1 and Qq ~ 0.2 
without a cosmological constant are shown in figure 4. The 
density profile (upper panel) is more concentrated in the 
open than in the flat case. This is similar to the behavior of 
density perturbations in the tophat model discussed above. 
The distribution and velocities of particles in the open model 
seem to be more evolved in time than in the flat model. 



where 



This relation works remarkably well for 0.1 < 57 < 1 and 
Si < 400. For generic configurations we assume that the re- 
lation ( |r^ ) is still valid. However we should take into account 
the fact that the "dimensionality" of the collapse affects the 
amplitude of the Q dependence; for example in the one di- 
mensional collapse, before shell crossing, the equations are 



4 n DEPENDENCE IN N-BODY 
SIMULATIONS 

We use N-body simulations to study the Q dependence under 
general initial conditions. These simulations are especially 
useful in orbit mixing regions where, according to (|l3|), the 
effect is most important. 

We ran two simulations having Ho = 1 and Clo — 
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Figure 4. The density profile (upper panel) and the "velocities" 
(lower panel) in the shell-crossing region in one dimension for 
Ho = 1 (solid) and Qq = 0.2 (dotted) with no cosmological con- 
stant. The initial density perturbation is a cosine wave symmetric 
about X = 0. 



Figure 5. The particle distribution in the low il (right) and f2 = 1 
(left) simulations. Slice thickness is lh~^Mpc. The lower panels 
focus on the group of "clusters" appearing near the centers of the 
upper panels. 



0.29 respectively. Both simulations were started from the 
same initial conditions. The initial conditions were gen- 
erated from the power spectrum for standard CDM with 
Ho = 50km s~^Mpc~^. Each simulation contained 128^ par- 
ticles in a cubic box of length 60h~^Mpc. The simulations 
were evolved until the linear rms density fluctuations in 
a sphere of 800km s~^ was 0.5. Both models have roughly 
the right small scale power as measured by the galaxy pair- 
wise velocity, but produce fewer rich clusters than observed. 
However, given the scale of our simulations, our choice of the 
power spectrum and the normalisation is appropriate for our 
purposes. A model with a higher normalisation would result 
in too much merging of smaller objects into a few larger ob- 
jects. Choosing a steeper power spectrum leads to a similar 
effect. 

The simulations were run using a modified parallel ver- 
sion (MacFarland et al. 1997) of Couchman's P'^'M code 
(Couchman et al. 1995) which uses explicit message pass- 
ing. The simulations had a softening parameter of 13.2% 
the mean particle separation and a mesh size of 512 in one 
dimension. They were run using 64 processors on the CRAY 
T3E supercomputer at the Computer Center of the Max 
Planck Society (RZG), Garching. 

Figure 5 shows the particle distribution in a slice of 
thickness Ih^^Mpc, in the two simulations. The left and 
right panels correspond to the flat and open model respec- 
tiveley. The lower panels zoom in on the "clusters" seen 
near the centers of the upper panels. On large scales (upper 



Figure 6. The two point correlation functions for the particle 
distribution in the two simulations. 



© 0000 RAS, MNRAS 000, 000-000 



6 Nusser & Colberg 



Figure 7. Densities in the open model vs. densities in the the 
flat model. The lower left panel shows densities after the CIC 
interpolation on a cubic grid of mean particle separation cell size. 
The other panels show CIC densities smoothed with a top-hat 
window of width Rs as indicated in the plot. 



panels) the two simulations arc remarkably similar. Some 
differences can be spotted in the lower panels. Clusters in 
the low fl simulation appear to be more concentrated and 
evolved. The differences between the the two simulations 
seem to be negligible on scales larger than lh~^Mpc. Indeed 
the mis value of the difference between the positions of the 
same particles in the two simulations is 0.25h~^Mpc and 
the largest difference is less than 1.5h~^Mpc. The correlar 
tion functions for the two simulations, plotted in figure 6, 
confirm the visual impression from figure 5. The correlation 
functions differ only on scales smaller than 1.6h^^Mpc. On 
"cluster" scales ^0.6h~'^Mpc the low f2 correlation func- 
tion is larger and, on scales 0.6 — 1.6h~^Mpc roughly cor- 
responding to the scale of infall regions around clusters in 
the simulations, it is smaller than the correlation function 
in the flat model. 

We now quantify the differences between the density 
and velocity fields. We first use the cloud-in-cell (CIC) inter- 
polation scheme to evaluate the density and velocity fields 
on a cubic grid of cell length equal to the moan particle 
separation. This produces a mass weighted average velocity 
on the grid points. We then further smooth the resultant 
density and velocity maps with a top-hat filter. In figure 7 
we plot the densities in the open vs those in the flat model 
for 4096 randomly chosen grid points. Even with only CIC 
smoothing on the grid scale, the correlation between the two 
density fields is very tight. For densities larger than 10 or 
so the densities in the open model are larger. The scatter 



Figure 8. The same as figure 6 but for peculiar velocities. The 
velocities in the open model were scaled by the factor /(f2o) = 
0.477. 



Table 1 . Moments of the density field in the two simulations after 
CIC and top-hat smoothing of width,_Rs, expressed in km/s. 







Ho 


= 1 


Oo = 


0.29 






o-<S 


S 




S 




CIC 


7.331 


6.222 


8.816 


6.374 


R 


s = 50 


6.313 


5.934 


7.397 


6.154 


Rs 


= 100 


3.331 


4.855 


3.603 


5.035 


Rs 


= 200 


1.693 


3.720 


1.752 


3.833 


Rs 


= 400 


0.878 


3.060 


0.891 


3.116 


Rs 


= 800 


0.440 


2.549 


0.443 


2.584 



almost vanishes when the density fields are smoothed with 
a top-hat window of width > 200km/s. Note however that 
for moderate densities (0 < 5 < 5), the densities in the flat 
model are slightly larger than the open model. This is not 
surprising since the general tendency is that matter flows 
out of regions with moderate densities into higher density 
regions. And since the open model is slightly more evolved, 
these moderate density regions are somewhat less dense in 
the open than in the flat model. In table 1 we list the values 
of the rms, as, and the reduced skewness, S =< 5^ > /aj 
of the density fields as a function of the smoothing scale. It 
appears, from the table, that the dependence of 5* on is 
stronger than what is predicted from second order perturba- 
tion theory (Bouchet et al. 1992). That theory is, however, 
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valid only for ct < 1. We can use table 1 to determine Afeff 
in ( ]l9| ) which relates ag in an open universe to that in an 
Q = 1 universe. A comparison of ( [l9| ) with table 1 suggests 
that Meff ~ 2. This value is reasonable since non-linear col- 
lapse configurations are likely to have pancake-like shapes. 
Recall that the simulations were stopped when the the lin- 
ear value of of as smoothed with a top-hat filter of width 
Rs = 800km/s was 0.5. The actual value computed from the 
simulation is very close to 0.44 in the two simulations. Thus, 
even though nonlinear effects are clearly important, the dif- 
ference between the as in the open and fiat simulations is 
negligible. We now consider the evolved velocity fields. Fig- 
ure 8 compares one of the components of the velocity fields 
in the two simulations. The velocity fields in the open model 
are scaled by the factor /(fio)- Even for large velocities and 
small smoothing widths, the velocity fields in the two simu- 
lations seem to be related by the factor /. A close inspection 
of the scatter plot for Rs — 400km/s reveals that the slope of 
the regression of v/ f on vi is slightly less than unity. This is 
because large velocities are generally associated with strong 
nonlinear effects which tend to spoil the scaling by /. 

We mentioned at the end of section 2, that the motion 
of particles in bound objects is independent of Q in terms 
of a time variable which corresponds to a velocity which 
is the peculiar comoving velocity divided by Q^'^ . Therefore 
we expect the velocity dispersion (rms velocity) in groups of 
particles identified in the simulation to scale, approximately, 
like QP'^ . To test this conjecture, we have used a friends- 
of-friends (f.o.f) algorithm (kindly supplied by A. Diaferio) 
to identify groups in the simulations. We then computed 
the one- dimensional velocity dispersion of particles in each 
group. Figure 9 (upper panel) shows the mean velocity dis- 
persion in groups as a function of the number of particles 
they contain. Velocities in the plot are scaled by Q,q'^ for 
the open model. It seems that this scaling works well. Given 
the uncertainties in identifying group members by the f.o.f 
algorithm, the deviations from this scaling for large groups 
are not significant. We conclude that while the (smoothed) 
average velocity of particles scales like /(fi), the rms ve- 
locity, roughly, scales like flP'^ . It is interesting to compare 
the abundance of groups in the two simulations. The lower 
panel of figure 9 is a plot of the abundance as a function of 
the number of particles for the two simulations. The abun- 
dance of groups is slightly higher in the open model. This is 
easy to understand, because groups in the open model are 
tighter than groups in the fiat model, the f.o.f algorithm nat- 
urally assigns more particles to them. Note that observations 
naturally provide the abundance of groups as a function of 
the mass. Since the mass of groups with the same number 
of particles is proportional to Qq, abundance curves, when 
plotted versus the mass, look significantly different in the 
two simulations. 



5 SUMMARY 

We have shown that gravitational dynamics of a pressure- 
less fiuid in an expanding universe is almost independent of 
the cosmological parameters. According to the equations of 
motion, expressed in terms of the linear growing mode, the 



Figure 9. Top: The mean 1-D velocity dispersion {rms velocity 
deviations) divided by Oj]'^ in groups vs. the number of particles 
they contain for the two simulations. The lengths of the error- 
bars give the 1-sigma scatter about the mean for the flat model. 
The scatter in the open model is similar. Bottom: The abundance 
of groups vs. the number of particles they contain. Groups are 
identified using a f.o.f algorithm. 

final structure in a low 51 model, with or without a cosmo- 
logical constant, is more evolved than in a fiat universe. We 
used toy-models and N-body simulations to investigate the 
effect of changing the cosmological background on the evo- 
lution of fiuctuations and, in particular, on the final velocity 
and density fields. The present density, when smoothed on 
scales slightly larger than cluster scale, is almost insensitive 
to the cosmological background. The background can affect 
the structure of bound objects (or halos). Halos, charac- 
terised by the same ratio of mass to background density, are 
more centrally concentrated in an open than in a fiat uni- 
verse. However, this effect is weak and is likely to depend on 
the initial power spectrum. On the other hand, the ampli- 
tude of the comoving peculiar velocities of particles strongly 
depends on fl. It is remarkable that the smoothed nonlinear 
velocity field scales with the growth factor, /, just as it does 
in linear theory. 

Since / depends very weakly on the cosmological con- 
stant, the final velocity field is mainly sensitive to il. There- 
fore the observed peculiar velocity and density fields in the 
nearby universe contain information only on Q (Lahav et al. 
1991). However, constraints on both 57 and A can, in prin- 
ciple, be obtained by measuring the clustering amplitude at 
different redshifts, for example via the correlation function 
(Lahav et al. 1991). As we have shown, it is rather difficult 
to detect signatures of the cosmological background in the 
structure of density fields. Fortunately, observations provide 
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the distribution of galaxies in redshift space. Thanks to the 
strong dependence of the velocity field on Q, the anisotropy 
of clustering in redshift space can provide a measure of f2. 
However, such estimates of fl involve an assumption on the 
relationship between the distribution of galaxies to that of 
the dark matter. Estimates of O independently of the galaxy 
distribution can be obtained from the observed peculiar ve- 
locity field alone (c.f. Dekel 1994). This makes peculiar ve- 
locity catalogs a very powerful too to constrain the cosmo- 
logical model. It is especially important for future peculiar 
velocity measurements to aim at larger sky-coverage and 
denser sampling rate. 
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